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Several prominent proposals have suggested that spins of localized electrons could serve as quan- 
tum computer qubits. The exchange interaction has been invoked as a means of implementing 
two qubit gates. In this paper, we analyze the strength and form of the exchange interaction un- 
der relevant conditions. We find that, when several spins are engaged in mutual interactions, the 
quantitative strengths or even qualitative forms of the interactions can change. It is shown that 
the changes can be dramatic within a Heitler-London model. Hund-Miilliken calculations are also 
presented, and support the qualititative conclusions from the Heitler-London model. The effects 
need to be considered in spin-based quantum computer designs, either as a source of gate error to 
be overcome or a new interaction to be exploited. 

PACS numbers: 03.67.Lx,75.10.Jm 



I. INTRODUCTION 

The exchange interaction between electrons has been 
studied since the early days of quantum mechanicsiS^ 
and has been reviewed in some classic references'*^ 
as well as textbooks^ Recently, a promising proposal 8 
has emerged to use the exchange interaction as a tun- 
able qubit-qubit interaction in a quantum computer, 
with the individual spins of electrons acting as qubits. 
To satisfy the conditions for constructing a universal 
quantum computer, the exchange interaction can ei- 
ther be supplemented with single-qubit operations^ or 
can be used by itself to construct a universal set of 
gates, in which case one encodes a logical qubit into 
the state of several spins AiSiiiiiS (This alternative to 
the standard universality scheme has been termed "en- 
coded universality"^) Motivated by the proposal of 
Loss and DiVincenzo^, there have been a number of 
studies of the one-particle and two-particle behavior of 
electrons localized on quantum dots within a quantum 
computer i 8 i 14 i 15 i 16 i 17 i 18 i 19 i 20 i 21 Here, we expand on our 
work 2 ^ considering the important situation of three or 
more coupled dots. We show how both quantitatively 
and qualitatively new effects can appear. These effects 
require consideration if one intends to make a quantum 
computer with more than two spins. 

The exchange interaction between two localized elec- 
trons arises as a result of their spatial behavior, but it 
can be expressed as an effective spin-spin interaction. In 
conditions of rotation symmetry (i.e. neglecting exter- 
nal magnetic fields, spin-orbit coupling, etc.), a purely 
isotropic form of this interaction arises, which is known 
as Heisenberg exchange: 

H cx = JS A • S B . (1) 

Here, S = (S x , S y , S z ) is a vector of spin- ^ angular mo- 



mentum operators, and A, B are indices referring to the 
location of each electron. (We take spin operators to be 
dimensionless in this paper - h is excluded from their 
definition.) This Hamiltonian has a spin-singlet eigen- 
state and degenerate spin-triplet eigenstatesi The quan- 
tity J is the exchange coupling constant, given by the 
energy splitting between the spin-singlet and spin-triplet 



states. 



J = e t 



(2) 



To date, studies of the exchange interaction in quantum 
computation have focused on the case of two quantum 
dots. 8 ! 14 ' 1 ^ 1 ?'^! 18 ! 1 ?! 2 ?] 21 Starting from the simplest case 
of two electrons in singly-occupied dots in the lowest or- 
bital state, systematic generalizations have been intro- 
duced and their effect on the exchange interaction stud- 
ied. In particular, researchers have analyzed the effect of 
double occupationji£ii2i2i higher orbital states^*^ and 
many-electron dotstf An accurate numerical study re- 
porting singlet-triplet crossing via magnetic field manip- 
ulation in a lateral double quantum dot can be found 
in Ref. |2jJ Neglecting spin-orbit coupling, these stud- 
ies have found increasingly accurate expressions for J, 
while focusing on the definition of Eq. (J2J . In the pres- 
ence of spin-orbit coupling both rotation and inversion 
symmetry are broken, and anisotropic corrections to H cx 
ariseM^ 

In this work, we undertake a study of the case of 
three or four electrons, each in a quantum dot. Once 
the system involves more than two electrons, simultane- 
ous multi-partite exchanges can occur. For three cou- 
pled dots containing three electrons, processes in which 
all three electrons exchange contribute to a quantitative 
correction to the value of J. We show explicitly that, 
for three identical dots arranged on the corners of an 
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equilateral triangle, the effective Hamiltonian can still 
be written using a Hciscnbcrg exchange interaction 

-ffspin = K + J (Sa ■ S B + S B • S c + S c • Sa) , 

but J is now found to be influenced by three-body ex- 
change matrix elements. For four coupled dots containing 
four electrons, the actual form of the interaction Eq. ||TJ 
changes due to four-body effects. For identical dots ar- 
ranged on the corners of a symmetric tetrahedron, the 
interaction takes the form 

H spin = K + J ]T S i -S j + J'[(S A -S B )(S C -S D ) 

A<i<j<D 

+ (S A • S C )(S B ■ S D ) + (S A ■ S D )(S B ■ S c )} 

where throughout this paper notation like A < i < D 
means that i takes the letter values A to D. Accord- 
ing to our Heitler-London (HL) calculations, the ratio 
\J'/J\ can reach 15% in in physically relevant parameter 
regimes. 

Four-body exchange terms have been discussed in 
other contexts - for example, in a perturbative treatment 
of the two and the three dimensional half-filled Hubbard 
models four-body interactions were shown to suppress 
the Neel temperature and the temperature of the para- 
ferromagnetic phase transition^LSI (see Appendix \Kfy . 
Here, we present a non-perturbative derivation of these 
terms, starting from a finite-dimensional Hamiltonian, 
and then highlight their significance for quantum compu- 
tation. Interaction Hamiltonian calculations like ours are 
of significance in various quantum computation contexts 
including (i) the encoded universality paradigm, where in 
the most efficient implementations several exchange in- 
teractions are turned on simultaneousl3«iii2*ii*i2iiii (quan- 
titative studies of parallel gate sequences-!- 1 in partic- 
ular may require revisiting in light of our results, as 
well as the "supercoherent qubits" method for reduc- 
ing decoherence^ where four and eight-spin interactions 
must be turned on simultaneously in order to enact quan- 
tum logic gates between encoded qubits); (ii) adiabatic 
quantum computing 29 , where the final Hamiltonian for 
any non-trivial calculation inevitably includes simulta- 
neous interactions between multiple qubits; (iii) fault- 
tolerant quantum error correction, where a higher degree 
of parallelism translates into a lower threshold for fault- 
tolerant quantum computation operations$22i2ii22i22i2i; 
(iv) the "one-way" quantum computer proposal^ where 
all nearest-neighbor interactions in a cluster of coupled 
spins are turned on simultaneously in order to prepare 
many-spins entangled state; (v) the search for physical 
systems with intrinsic, topological fault tolerance, where 
systems with four-body interactions have recently been 
identified as having the sought-after properties^ 

We begin with a general description of a finite- 
dimensional effective spin Hamiltonian in Section [H] 
(This is compared to the standard, perturbative deriva- 
tion in Appendix lAll Section IlIII shows how to compute 
the parameters in the effective spin Hamiltonian, with de- 
tailed consideration of the two-electron, three-electron, 



and four-electron cases. We introduce a specific model 
and calculate the parameters quantitatively, for three and 
four electrons, in Section llVl Appendices [B] and El con- 
tain relevant technical details. 



II. ELECTRON-SPIN-OPERATOR 
HAMILTONIAN 

In this section we present general arguments concern- 
ing the form of the effective spin Hamiltonian, as it arises 
from n localized electrons interacting via the Coulomb 
force. We start with the familiar electronic Hamiltonian 

n i 2 

i=l i<j 12 - 71 

n 

= ^2h(Ti) +^2<w(r i ,T j ), (3) 

i=i %<j 

where the first term is the kinetic energy, the second is 
the confining potential, and the third is the Coulomb in- 
teraction. The confining potential V(r) contains n energy 
minima, which give rise to the n dots. To understand the 
dynamics of n electron-spin qubits in n quantum dots, it 
is desirable to eliminate the spatial degrees of freedom, 
leaving an effective Hamiltonian composed of electron- 
spin operators only. 

The first step in changing Hamiltonian © to an 
electron-spin-operator Hamiltonian is to fix a basis. We 
first consider the case of two electrons in two dots, labeled 
A and B. We do not allow for double occupancy of a dot 
and consider only a single low-energy orbital per dot la- 
beled as 4>a(y) = ( r l^4) and <^s(r) = (r\B) . Electrons 1 
and 2 occupy these low energy orbitals. Each electron 
can have spin-up or spin-down, hence each electron can 
represent a qubit. A state with, for example, electron 1 in 
orbital B with spin-up and electron 2 in orbital A with 
spin-down is represented as \BA) Since electrons 

are fermions, this state needs to be antisymmetrized; the 
full state of the two electrons takes the form of a Slater 
determinant 

|*UT)) = (\AB) |4T> - \BA) |n)) cx a A[ a B] |vac) . (4) 

Note that the order of spins in the state label ^(IT) in- 
dicates that the electron in orbital A has spin down and 
the electron in orbital B has spin up. In Eq. we 
introduced second-quantized notation (ignoring normal- 
ization), with a A ^ creating an electron with spin down 
in orbital A and a g ^ creating an electron with spin up 
in orbital B. The four states \9(sa,sb)) form the two- 
electron basis. The same procedure applies to three elec- 
trons in three dots. There is again a single low-energy 
orbital per dot, labeled as </u(r) = (y\A), ^s(r) = (v\B), 
and </>c( r ) = ( r |C) for dots A, B, and C respectively. 
Electrons 1, 2, and 3 occupy these low energy orbitals. A 
state with, for example, electron 1 in orbital B with spin- 
up, electron 2 in orbital A with spin-down, and electron 3 
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in orbital C with spin-up is represented as \BAC) |Tlt) • 
This state is then antisymmetrized so that the full state 
of the three electrons is the Slater determinant 

|*UTT)) = (\ABC) UTT> + \BCA) |TU) + \CAB) |UT) 
- - \BAC) |UT) - \ACB) HIT) - \CBA) |TU)) 
oc a^a^a^ |vac) . (5) 

The order of spins in the state label ^(ITT) indicates 
that the electron in orbital A has spin down, the electron 
in orbital B has spin up, and the electron in orbital C 
has spin up. The eight states \^>(sa,sb,sc)) form the 
three-electron basis. 

The general case of n electrons gives 2™ fully antisym- 
metrized basis vectors of the form 

|tf (sa, sb, ... , s z )) = SpP[\AB ...) \ SA s B ...)], (6) 
p 

where the sum runs over all permutations P of both or- 
bitals and spins, and 8p = 1 (—1) if the permutation is 
even (odd). In this basis, the Hamiltonian takes the 
form of a 2" x 2" Hermitian matrix. Like any 2" x 2™ 
Hermitian matrix, the Hamiltonian can be written as a 
sum 

3 

H spin = ^(A) (g> (TjiB) (g> ■ ■ ■ (7) 

i,i,~=o 

of Hermitian spin matrices of the form <Ji(A)(8)<Tj(B)®. ■ ■ 
each multiplied by a real coefficient hj,.-.. Here, <Ji(p) 
denotes the Pauli matrix <jj acting on the electron in 
dot p, with i = 0, 1, 2, 3 and with ag equal to the iden- 
tity matrix. There are n factors in the tensor product 
<Ji{A) (8>CTj (B) ® • • • , so that it can be written as a 2" x 2" 
matrix, and there are similarly n subscripts on the coef- 
ficient This decomposition J7J) into spin matrices 
produces an effective electron-spin Hamiltonian that con- 
veniently describes the dynamics of n qubits. 

The procedure we have just described is framed within 
the Heiter-London approximation^ The approximation 
consists in neglecting excited states and has been criti- 
cized on the grounds that it does not produce the correct 
asymptotic behavior in the limit of very large distances^ 
However, in the context of our system of interest, this 
asymptotic limit is not a concern, and moreover, recent 
studies have verified the utility of the approximation in 
the case of large (but not infinite) inter-dot separation^ 
We will thus proceed with the HL approximation, which 
has the advantage of conceptual simplicity and physical 
clarity. In the three electron case, we show that Hund- 
Miilliken (HM) calculations, in which double occupation 
is permitted, support the conclusions of our HL results. 

Symmetry considerations fundamentally constrain the 
form of the electron-spin Hamiltonian. The coordinate 
system used to define f and J, is arbitrary if there is no 
spin-orbit coupling and no external magnetic field. In 
this case, the effective spin operator Hamiltonian has 



rotation, inversion, and exchange symmetry. The co- 
efficients hj,--- in J7J) are strongly constrained by this 
symmetry. The Hamiltonian can only be a function of 
the total spin squared = (Sa + Sb + . . . ) 2 , where 
S A = ±{a 1 (A)x + a 2 (A)y + <T 3 (A)z)®ao{B)(E><7 () (C)(E)- ■ • , 
S B = ao{A)^^{a 1 (B)x+a 2 (B)y+a 3 (B)z)^a a (C)^- ■ ■ , 
etc. A scalar such as Sa ■ (Sb x Sc) cannot appear in the 
Hamiltonian because of inversion symmetry. We must 
have 

-ffspin = Lq + L\S T + L 2 (S T ) 2 + . . . (8) 

where Lq, L%, L 2 , ■ ■ ■ are real constants with dimensions 
of energy. The constant Lq is an energy shift. The term 
proportional to L\ gives rise to the familiar Heisenberg 
interaction. Here we see that in principle higher order in- 
teractions may be present in the spin Hamiltonian, start- 
ing with a fourth order term proportional to L 2 . In this 
highly symmetric situation, the eigenstates of the spin 
Hamiltonian are clearly just eigenstates of S^- 

III. COMPUTATION OF THE SPIN 
HAMILTONIAN PARAMETERS 

To compute the values of Lq, Li, L2, ■ ■ . we consider an 
eigenstate |^) of S 2 ^, with known eigenvalue St(St + 1)- 
If there are n electrons in the system, we write \^) = 
)• To proceed, one (i) computes the expectation 
value of the effective spin Hamiltonian (JSJ in this state, 
(ii) computes the expectation value of the spatial Hamil- 
tonian © in this state, and then (iii) equates the two 
expectation values: 

(*iff spin i*> - (mm- (9) 

This procedure is repeated for all eigenvalues of S 2 ^, thus 
generating a set of linear equations for the parameters 
Lq, Li, Li, . . . , in terms of matrix elements of H between 
different orbital states. For n electrons the number of 
distinct eigenvalues of S\ is |_§ J + 1 (where [§J denotes 
the greatest integer less than §), so this is the maximum 
number of distinct energy eigenvalues of the Hamiltonian 
©■ Thus, the coefficients L m for < m < [f J + 1 have 
enough degrees of freedom to completely and uniquely 
specify the matrix @ ; without loss of generality, we can 
set L m = for to > [f J + 1- We are led to [f J + I 
coupled linear equations for the non-zero L m parameters. 
In the case that n is even St takes on the integer values 
0, 1, . . . , n/2. In the case that n is odd St takes on the 
half-integer values 1/2, 3/2, . . . , n/2. Then 

LfJ 

(n T \H sP in\n T ) = j2L m (s T (s T +i)r m 

Having completed step (i) of our program, we now turn 
to step (ii), the calculation of (*g T |ff|*g r ). We make 
this calculation separately for the cases of two, three and 
four electrons. 
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A. Two Electron Case 

As a simple illustration of our procedure we rederive 
the well-known result for two electrons: the exchange 
constant equals the difference between the (degenerate) 
triplet states and the singlet state. The spin singlet 
(St = 0) and spin triplet (St = 1) states have eigen- 
values of S T equal to and 2 respectively. Thus, the 
Hamiltonian JSJl can only have two distinct eigenvalues, 
and we need to solve |_§J +1 = 2 equations for Lq and L\. 
A convenient St = 1 eigenstate is the normalized state 
1*3=^) = JV|tf(TT)) = M(\AB)\]]) - \BA)\\V) = 
J\f(\AB) - \BA)) |tt) • The normalization constant J\f has 



which can be evaluated once a choice of orbital states 
is specified; we do this in Section \W\ below. To com- 
pare with equation we note that this equation spec- 
ifies the triplet energy e t = Lq + 2L\. A second 



giving the singlet energy e s = Lq. To exhibit the ex- 
change coupling explicitly, we rewrite the Hamiltonian 

as 

H S pin = Lq + LiS\ + LiS^ + 2LiSa ■ Sb 

= K + JS A -S B . (13) 

where K = L + (3/2)L 1 and J = 2Li. Expression @ 
follows when we note that L\ = i (e t — e s ). 

B. Three Electron Case 

1. Heitler London Model 

In the three electron case, the possible values that 
the total spin can take are St = 1/2 (with two, two- 
dimensional eigenspaces) or St = 3/2 (with a four- 
dimensional eigenspace). We therefore again need to 
solve L|J +1 = 2 equations, and it is sufficient to keep 
only two constants Lq and L\ in H spin , setting Li and the 
rest to zero. As a convenient state with known St = 3/2 
we take the normalized state | ^ | / 2 ) ^ I^CTTT))) so that 



the value 

N = ({AB\AB)+(BA \BA)-(AB \BA)-(BA \AB))~ 1/2 . 
Inserting this state into Eq. © yields 

The spin Hamiltonian's expectation value is immediately 
found to be Lq + 2Li, as can be seen from Eq. Hl()|l . Ex- 
panding out the spatial Hamiltonian's expectation value 
gives 



(11) 

I 

equation is found from the St = state ^^ig) = 

M(MU)) ~ \mW) = AT(I^B) + \BA)) (|U) - ut», 

which leads to 



(12) 

I 

the energy is E 3/2 = (^l /2 \H svin \^\ /2 ) = L + 
We use \^ 3 1/2 ) cx ^ (|*(UT)) - l*(itT)» a « a normal- 
ized state with known St — 1/2 , for which the en- 
ergy is E 1/2 = (*? /2 |ir.pi»|*? /2 ) = Lq + Li(£)(§). 
Then equating expectation values of Hamiltonian ({3*)l 
and Hamiltonian ©, i.e., requiring (^g -ffspin ) = 
(*| I H ) for each of our states St = 1/2 and 
St — 3/2 as in Eq. we can solve for Lq and L\. 
To do so we need to obtain more explicit expressions for 
(*| T | H |*| T )- We assume that (f> A ( r), B (r) , and 
4> c (h are real and satisfy (A\A) = (B\B) = (C\C) and 
(A\B) = (A\C) = (B\C) (this is consistent with our orig- 
inal assumption of rotation, inversion, and exchange in- 
variance). First, let us normalize I ^3/2) : 

|*! /2 ) = N[\ABC) + \CAB) + \BCA) 

-\BAC) - \CBA) - \ACB)\ |fTT) , 

where the normalization constant Af is given by 



>/6 (p 3 + 2po - 3pi) 



Lq 



2Li = 



(ab \ h \ab) + (ba\ h \ba) - (ab\ h \ba) - (ba\ h \ab) 
Jab \ab) + (ba \ba) - (ab \ba) - (ba \ab) 



Lq = 



(AB \ H \AB) + (BA\ H \BA) + (AB\ H \BA) + (BA\ H \AB) 
{AB \AB) + {BA \BA) + {AB \BA) + (BA \AB) 
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c 

o 



These equations give Lq and L\ in terms of the and 

To compute the usual exchange coupling, it is useful 
to rewrite H sp i n as 



where 



A<i<C A<i<j<C 

= K + J(S A -S B + S A -S C + S B -S C ). (16) 



K 



Lr, 



-M 



J = 2L X . 



(17) 
(18) 



Solving for the exchange constant J — 2L\ we find finally: 



FIG. 1: Matrix elements relevant to three electron case. Ar- 
rows indicate transition from localized state on initial dot to 
localized state on final dot. 



The quantities p 3 , pi, and po are given by 

p 3 = (ABC\ABC), 

which is an overlap integral when all three electrons re- 
tain the same state in the bra and kct, 

pi = (BAC\ ABC) = (CBA\ ABC) = (ACB\ ABC), 

which is an overlap integral when one electron has the 
same state in the bra and ket, and 

Po = (CAB\ ABC) = (BCA\ ABC), 

which is an overlap integral when zero electrons have 
the same state in the bra and ket - all three electrons 
change their states. In evaluating the matrix element 
(^g, 2 \H\^>l, n ) we use the notation 



"3/2/ 

e = (CAB\ H \ABC) = (BCA\ H \ABC) 
ei = (BAC\H\ABC) = (CBA\H \ABC) 

= (ACB\H\ABC) 
e 3 = (ABC\H\ABC), 



where the physical interpretation is that Ck involves 3 — k 
electrons exchanging orbitals (Fig. [IJ. 

Computing the expectation value of H in the state 
/ 2 ) then leads to the result: 

E 3/2 = L + -L 1 ^ e 3 + 26 -3 £l 
3/2 ° 4 1 p 3 + 2p Q ~3 Pl V ' 

For the case St — 1/2, using / 2 ) an analogous cal- 
culation yields: 



p - T + 3 t - 63 ~ e ° 
4 p 3 - po 



(15) 



J — g {E3/2 - E1/2) ■ 

A couple of comments are in order concerning this result. 
First, the energies E 3 / 2 ,E 1 / 2 can be calculated once the 
orbitals are specified, as we do in Section Hvl below. We 
see that, similar to the two-electron case, the physical 
interpretation of the exchange constant is that (up to a 
multiplicative factor) it is given by the energy difference 
between the St = 3/2 and St = 1/2 states. Second, note 
from Eqs. I|14f) and (|15f) that the value of the exchange 
constant J is determined in part by the "three-electron- 
exchange" terms of the form po — (CAB \ABC) and 
eo = (CAB\H\ABC). It is apparent that such terms 
involve a cooperative effect between all three electrons 
and hence cannot be seen in two-electron calculations. 
It follows that the presence of the third electron quanti- 
tatively changes the exchange coupling between the other 
two electrons. 



2. Hund-Miilliken Model 

We have have been working within the HL approxi- 
mation in which there is one orbital per quantum dot 
occupied by a single electron. To check its physical va- 
lidity, we make three-electron computations within the 
HM approximation as well, in which double occupation 
of quantum dots is permitted. This leads to a total of 
8+12 = 20 basis states in the three spin case (2 3 = 8 
from the HL basis and 3x2x2 = 12 double-occupation 
states) . In the HL approximation, the 8 states divide into 
a degenerate four-dimensional S = 3/2 subspace with 
energy £3/2 and a degenerate four-dimensional S =1/2 
subspace with energy £1/2. In the HM case, the degen- 
erate four-dimensional 5 = 3/2 subspace is unaffected 
by the new double-occupation states which must all have 
S = 1/2; the energy of these four S = 3/2 states remains 
£ 3 / 2 . (The S = 3/2, S z = 3/2 state has three spin-up 
electrons and so the Hamiltonian cannot mix it with any 
other state. Since the other S = 3/2 states are related 
by a rotation of the arbitrary spin axis, they must be 
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eigenstates of the Hamiltonian with the same energy.) 
The 12 double-occupation states enlarge the S = 1/2 
subspace, which becomes 16-dimensional and has a non- 
trivial spectrum. 

In the HM case, the decomposition JJJ is no longer 
meaningful since the basis states do not necessarily have 
one spin per quantum dot. This complicates the com- 
putation of the eigenspectrum of this 16 -dimensional 



space. First, we note that the projection S z (the number 
of spin-up electrons) is still a good quantum number since 
the Hamiltonian J3J cannot mix two states with different 
numbers of spin-up electrons. The 16-dimensional sub- 
space therefore splits into two degenerate 8-dimensional 
S z = ±1/2 subspaces. The S z = 1/2 subspace consists 
of two HL states and six double-occupation states anal- 
ogous to (JSJ: 



J 



1*5 



ex -L(|*(UT)) 
,t J J 



l*UTT)»,l*? 



/2 > oc — |*(TU)) 



l 



V6 



(i*(UT)> + i*an)», 



|*aab(TIT)> cx a\ ]a \^ B] I vac) , |*aac(UT)) , |*bba(UT)) , |*bbc(UT)) , |*cca(UT)) , |*ccb(UT)> ■ (19) 

I 



One can construct the 8x8 Hamiltonian in this sub- 
space and diagonalize it. The eigenstates exhibit degen- 
eracies arising from the symmetry of the Hamiltonian 
under the exchange of a pair of dots. Assuming that 
our dots are all equivalent, there are three dot-pair ex- 
change operators that commute with the Hamiltonian: 
Ea,b that exchanges dots A,B, Eb,c that exchanges 
dots B, C, and Ec,a that exchanges dots C, A. For in- 



stance, E A ^ B |*aac(T4T)> = |*bbc(UT)>- We can re- 
quire that the eigenstates of the Hamiltonian also be 
eigenstates of Ea,b ° r Eb,c ° r Ec,A- Using our 8 states 
(|19fl . it is possible to construct 2 linearly independent 
states that are simultaneous eigenstates of all three ex- 
change operators. The two (unnormalized) eigenstates 
are 



J 



(|*WUT)) + l*AAc(Ut)) + |*BB.4(TIT)> 

(I^ b (tit)) - i*aac(T4T)) - \*BB A (nm 



I*BBC(TIT)> 

I*sbc(T4T)) 
I 



|*cca(UT)) + |*ccb(UT)» and 

|*cca(UT)> - |*ccu(UT)»> 



with eigenvalue +1 and — 1, respectively. Each such state 
turns out to be an eigenstate of the Hamiltonian with its 
own non-degenerate energy. The remaining 6 members of 
the 8-dimensional subspace are not simultaneous eigen- 
states of all three exchange operators Ea,b and Eb.c 
and Ec,a- To ensure that we can nevertheless choose 
the eigenstates of the Hamiltonian to be simultaneous 
eigenstates of Ea,b or Eb,c or Ec,Ai the energy eigen- 
states occur in degenerate pairs that can be superposed 
as desired to form eigenstates of the exchange opera- 
tors. When the parameters of the spatial Hamiltonian 
@ make double occupation energetically expensive, one 
of the degenerate pairs will be low in energy and will 
consist mainly of the HL states |^f/ 9 ) and |3^/ 2 )- In 
this way, the HM calculation reduces to the HL result 
plus high energy double-occupation states, and 116fl still 
describes the low-energy spin dynamics. 



3. Unequal Coupling 

We emphasize that Eq. I|16H was derived assuming rota- 
tion, inversion, and exchange symmetry. Exchange sym- 
metry, in particular, is broken whenever there is unequal 
coupling between dots, and then the Hamiltonian can 
involve more constants. This situation is realized when 
the dots are not all equidistant, or when they have been 
shifted electrically, as in the case of dots defined by elec- 
trodes creating confinement potentials, 38 or when there 
are unequal tunneling barriers between different dotsJ^ 
For instance, in the case of three unequally coupled dots 
the Hamiltonian will have the form 



ffspin - K 



JabSj 



Jbc&2 ■ Sc + Jac^a ■ Sc 



if we still assume rotation and inversion invariance. (An 
external magnetic field, which has been shown to be in- 
strumental in changing the sign of J in the case of two 
dotSfi^±& could lead to a very different Hamiltonian. It 
would break rotation symmetry, introducing operators 
into -ffspin like S^.) This three electron Hamiltonian 
commutes with the z-component of the total spin op- 
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erator St, so they can be simultaneously diagonalized. 
The state |^(fj.l)) has S T = —3/2 and energy eigen- 
value K + {Jab + Jbc + Jac)/^- It is found to be 
degenerate with ^(|*(UT)> + |*(IU)) + |*(TU)»> a 
state with S T = —1/2. There are two remaining S T = 
— 1/2 eigenvectors, which have the (unnormalizcd) forms 

( J 5°:^ J l*(ut)) + J ^:i^ J > (m)) + i*(tu)»> 



Jab— J ao 

l+ll 



( J %o-j£ J l*(UT)) +3g§^^l*aTi)) + l*(TU)» 

and have energies a + J, a — J respectively, where a 
-3(J AB + Jbc + Jac)/2 and J = (J^ + J| c 



Jic 



JacJbc ~ JabJac - JabJbc) 1/2 - The remaining four 
energy eigenvectors, with Sf. = 3/2 and Sj, = 1/2, can 
be obtained from these four by inversion. From these re- 
sults it is possible to derive equations analogous to 114|) - 
ff 1 81) in the case when Jab, Jbc, and Jac are n °t equal. 



C. Four Electron Case 

In the case of four electrons, the effective Hamiltonian 
again takes the form (fHJ . Since four electrons can have 
St = 0, St = 1, or St = 2, we must keep three con- 
stants Lo, L\, and Li in H sp i n . It follows immediately 
that -ffspm includes terms of the form ^(S^-SbXSc-Sb) 
and permutations. Unless Li happens to vanish, the pres- 
ence of a fourth electron introduces a qualitatively new 
4-body interaction as well as a quantitative change in the 
exchange coupling between the other electrons. 

We now calculate Lq, L\, and Li just as we calculated 
Lq and L\ for three particles. Let us define 



Po 


= (BADC\ ABCD) 


Q) 


= (BADC\ H\ABCD) 


p'o 


= (DABC\ ABCD) 




= (DABC\ H \ABCD) 


Pi 


= (ADBC\ ABCD) 


ei 


= (ADBC\ H\ABCD) 


P2 


= (BACD\ ABCD) 


£2 


= (BACD\ H\ABCD) 


Pi 


= (ABC D\ ABCD) 


C4 


= (ABCD\H\ ABCD) 



where the subscript indicates how many electrons retain 
the same state in the bra and the ket, just as in the 
three electron case. The terms eo and e' involve four- 
body effects: eo involves two pairs of electrons exchang- 
ing orbitals, and e' involves all four electrons exchanging 
orbitals cyclically (Fig. |2J • 

A convenient state to use for St — is 

^(i*(uu)> - i*(tut)) 7 i*am)> + i*atit)», 

keeping in mind the definition ©. After normalization, 
this state yields the singlet energy 



En — Lq — 



e 4 - 4ei + 3e 
P4, - 4pi + 3p 



A convenient state to use for St = 1 is = 

JV(|*(UU)> + l*(TUT)> - l*(ltti)> - l*UUT)». 

This state, after normalization, yields the triplet energy 
E 1 =L + 2 Ll + 4L 2 = £4 '„ 2£2 ' £0 + 2 > / (21) 



Pi - 2p 2 - po + 2p' 



© 




p o' e o 




® 


© 




FIG. 2: Selected matrix elements relevant to four electron 
case. Arrows indicate transition from localized state on initial 
dot to localized state on final dot. 



Finally, a convenient state to use for St — 2 is ^(TTTT)) 
We find for the quintet energy 



£4 - 6e 2 + 8ei + 3e - 6e 

Pi - GP2 + 8pi + 3p Q - 6p' ' 

(22) 



E 2 = Lq + 6Li + 36L 2 = 

Solving, we have 

Lq = E 

L l = ~( i E 2 -9E 1 +8E ) 
Li = ^(E 2 -3E 1 +2E ). 

We would like to exhibit interaction constants explic- 
itly in the spin Hamiltonian. We have St = Y^f=A Sii so 
that 

S rp — — 3 1 ~\~ 2 ^ ^ S i ' S j 

i<j 



(20) 



while 

(S T ) 2 = (3/+ £ S'-^) 2 

A<i^j<D 

and it can be shown that 
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i^J k^l 

= (S i -S i )(S fc -S i ) + 4 £ (S i -S J -)(S,--S0+2^(S i -S i )(S i -S J -) 

2 (S i -S J )(S fe -S0 + 4[i^Si-S i ) + 2(~-~$> i -S i ] . 



We are led to 
27 

(S|) 2 = y + 75>.S,+ ^ (Si-S^CSfc-S,). 

The spin Hamiltonian can now be written as 
# S pi„ = K + JJ2 Si ■ Sj + J'[(S A • S B ) (S c • S D ) 
+ (S A ■ S c ) (9 B ■ Sd) + (Sa ■ S D ) (S b • S c )]. 



where 



K = L + 3Li 



27L 2 



-2£ + 9E-L + 5E 2 



J = 2L 1 + UL 2 = 



2 16 
-2E - 3Ei + 10E 2 
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J' = 8L 2 = 



—2Eq — 3Ei + E 2 



(23) 
(24) 
(25) 



Generically, J' does not vanish, and four-body interac- 
tions arise. The physical interpretation of the exchange 
constants as simple energy differences between different 
spin multiplets is now lost: we find energy differences 
with numerical coefficients that are not intuitively obvi- 
ous. 

Of central physical importance to us is the relative 
sizes of the coefficients J and J'. This is studied in the 
next section, where a HL calculation suggests that J' 
is substantial in comparison to J in physically impor- 
tant regions of parameter space. We also find that both 
coefficients are affected by three (pi, ei) and four-body 
exchanges (p , p' 0) e , e' ). 

In the general case of 2n electrons, two-body, four- 
body,..., 2n -body interaction terms appear in the Hamil- 
tonian. Computing the strengths of the interactions for 



larger n is a topic of interest, but we do not address it 
here. One expects the strengths of the terms to decrease 
with the number of bodies involved. 



IV. MODEL POTENTIAL CALCULATIONS 

To compute the values of the L i: we select the following 
specific form for the one-body potential in © : 



V(t) = 



1 



A| 2 |r- 



-B| 2 |r-C| 2 |r-D| 



(26) 



2(20 6 

This potential has a quadratic minimum at each of 
the vertices of an equilateral tetrahedron A = (0,0,0), 

B = (2^,0,-2/^1), C = (-l^Jl,l,-2lJf), and 

D = {—l-\f\, —I, — The distance between vertices 
is 21. We select a potential with four minima so that 
it can be used in the four electron case without modi- 
fication. This facilitates comparison between the two-, 
three-, and four-electron cases, and the extra minima do 
not influence the two- and three- electron cases in any 
significant way. 

At vertex A, we define the localized Gaussian state 

/ / \ , | A \ /mw \ 3 / 4 / mu) 
Mr) = (r\A) = (— J eX p(- — |r-A| 

which is the ground state of the quadratic minimum at 
that vertex. We define localized states similarly for the 
other vertices. 

The following one-body Hamiltonian matrix elements 
are needed to evaluate the coupling constants in H sp i n 
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(A\A) = 1 (27) 



(A\h\A) 


P 2 


^m^v - A\ 2 \A) + (A\V(v) -^mu) 


o|r-A| 2 |A) 




= Tiu) 


"3 

2 4 


-^-(63x h - 3 + 280xr 2 + 320a; 6 - 1 ) 
2048 




(A\B) 


= e' Xb 








(A\h\B) 


= Tiojo 


"3 

r 


^ + „ n , (945xr 3 + 1680a;r 2 + OSe^ 1 1216 880x b )e- Xi 
2048 



In these equations, we have added and subtracted a har- 
monic oscillator potential from the one-body Hamilto- 
nian h for case of calculation. The dimcnsionless tunnel- 
ing parameter Xb is the square of the ratio of the interdot 
distance 21 to the characteristic harmonic oscillator width 



2y/h/(mu) ): 



Xb — mu) l 2 /h, 



which is also the ratio of the tunneling energy barrier 
^muj 2 l 2 to the harmonic oscillator ground state energy 



by 



Matrix elements of the Coulomb interaction are given 



(FG\w\UV) = fiujo 
(FG\w\UV) = hiu 



2x c \J Xb 

4 



2\/2 



-X c yJXb& 



f + u - g - V 

(l/4)(|f-u| 2 + |g-v| 2 ) 



-(l/4)(|f-u| 2 + |g-v| 2 ) erf 



|f + u-g- v| 



If + U- 



2\/2 



, |f + u-g-v|^ 



v| = 0. 



(29) 



In these equations, the dimensionless parameter x c is the 
ratio of the Coulomb energy e 2 /(n2l) to the harmonic 
oscillator ground state energy -^hu> : 

x c = e 2 / (nlhu) ). 

The symbols F, G, U, and V take values from the 
set {A, B,C, D}. The lower case vectors are defined 

by f = V^F F = ^ F > etc - The symbol erf(x) = 
-7= fn e~ s ds denotes the error function. 



A. Two electrons 

In the case of two electrons, we assume that two of 
the potential minima of l|26|) are occupied - there is 
an electron at A = (0,0,0) and an electron at B = 

(2Zy^, 0, — 2/^/|). In order to compute Lq and L\ 

from Eqs. and H12[> above, we require only the ma- 
trix elements (AB\ H \AB) = 2{A\h\A) + (AB\w\AB), 
(AB\AB) = 1, (AB\H\BA) = 2(A\h\B) (A\B) + 
(AB\w\BA), and (AB\BA) = (A\B) 2 . (We have sim- 
plified using the fact that (^4|/i|^4) = (B\h\B) and using 



the fact that the wavefunctions are real.) Once Lq and 
L\ have been computed, it is straightforward to obtain 
K = L Q + (3/2)Li and J = 2L V 

A plot of the energy shift K as a function of x b (the 
tunneling energy) and x c (the Coulomb energy) is shown 
in Fig. |3| in units of hui . Following Ref. |8j, we estimate 
realistic values for x b and x c by considering the case of 
GaAs heterostructure single dots. An estimated value for 
Xb is Xb = mio 2 l 2 1 '(hu) ) ~ 1, since the harmonic oscil- 
lator width 2y / h/ (miOo) should be approximately equal 
to the distance between dots 21 in a quantum computer. 
The parameter x c = e 2 / (nlhio ) m 1.5 taking k = 13.1, 
hu> = 3 mcV and i(, w 1, Note that the energy K in- 
creases when the one-electron tunneling barrier energy 
decreases and the Coulomb interaction energy increases 
(i.e. for small Xb and large x c ). 

In Fig. we plot the exchange-interaction constant J 
as a function of x b and x c - The plot generally indicates 
that J increases as the tunneling barrier decreases (xb 
smaller), an intuitively reasonable result. However, when 
the Coulomb interaction is particularly strong, as when 
x c —>■ 5, this trend is violated: J develops a negative 
minimum at x& ~ 1.5. The reason is that the St = 3/2 
state has a totally antisymmetric spatial wavefunction, 
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FIG. 3: Plot of K as a function of dimensionless tunneling 
barrier Xb and Coulomb energy x c in the case of two interact- 
ing electrons. 




FIG. 4: Plot of J as a function of xt and x c in the case of two 
interacting electrons. When J < the HL approximation has 
broken down. 



while the St = 1/2 state does not. The antisymmetry 
tends to reduce the Coulomb repulsion energy between 
electrons while increasing the one-electron-tunneling en- 
ergy. When parameters are tuned to make the Coulomb 
repulsion important, the energy of the St — 3/2 state 
dips down, eventually decreasing below the energy of the 
St = 1/2 state. This leads to J < 0. The negative 
value of J signals the breakdown of the HL approxima- 
tion in this region. The exact two electron ground state 
is known^S to have St = 1/2, while J < would imply 
an St = 3/2 ground state. The HL representation of the 
St = 1/2 state is simply too rigid to represent the exact 
ground state when interactions are strong. The inflex- 
ibility of the HL wavefunctions should be kept in mind 
when there are strong interactions in the three electron 




FIG. 5: Plot of X as a function of xt and x c in the case of 
three mutually interacting electrons. 

and four electron CBjSC, c.IS well. Fortunately, in our region 
of greatest interest, Xb ~ 1, x c sa 1.5, HL results should 
be meaningful. Even then, however, it should be kept in 
mind that the barrier between minima of the potential 
(|26|l is shallow, and so the calculation will become in- 
creasingly inaccurate as the minima get close together. 
Alter native nume rical methods can be found, e.g., in 
RBfs. ll4ll5llril7t 

B. Three electrons 

In the case of three electrons we assume three of the po- 
tential minima in f l26fl are occupied at, say, A = (0, 0, 0), 

B = {2lJ\,Q,-2lJfj, and C = {-lJ\,l,-2lJfj. 

The electrons are therefore arranged at the corners of an 
equilateral triangle, and the minimum at D is unoccu- 
pied. To solve 1|14[> and p5|l . we need to evaluate three- 
body matrix elements. Details are given in Appendix IBI 

A plot of the energy shift K as a. function of Xb (the 
tunneling energy) and x c (the Coulomb energy) is shown 
in Fig. [3 in units of hu> . The plot's shape is quite 
similar to that of Fig. |3J Fig. displays the change 
AK given by subtracting from K the value that K 
would take if the three-electron swap matrix elements 
e = {CAB\H\ABC) and p = (CAB\ABC) were zero. 
The axis directions are reversed in this plot to make its 
shape easier to inspect. The figure shows that AK is 
most important when one-electron tunneling barrier en- 
ergy the Coulomb interaction energy are small in magni- 
tude (small Xb and small x c ). 

In Fig. [7| we plot the exchange-interaction constant 
J as a function of Xb and x c . As in the two-body case, 
the HL approximation becomes suspect in the regime of 
strong interactions, so the region of negative J could just 
indicate its breakdown. However, the exact ground state 
for three or more electrons is not as well understood as 
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FIG. 6: Plot of AK as a function of xt and x c in the case 
of three mutually interacting electrons. Axis directions are 
reversed from the preceding figure. 




FIG. 7: Plot of J as a function of Xb and x c in the case of 
three mutually interacting electrons. 



that for two electrons^ the result J < cannot be ruled 
out a priori. 

Fig.[8]shows the change A J given by subtracting from 
J the value that J would take if the three-electron swap 
matrix elements eo and pg were zero (note that the axis 
directions are flipped to make the plot clearer). Com- 
paring the scales of Figs. [7| and 00 one finds that the 
three-electron swap matrix elements can have a powerful 
influence on J. 

To complement our HL results, we have computed the 
HM spectrum. For reasonable parameter values (xb — 
1.0, x c — 1.5), we have found that the lowest four states 
of the 16-dimensional 5=1/2 subspace are degenerate 
and have an energy (that we can call -Ei/2,hund) that is 
well-separated from that of the remaining 12 states with 
5 = 1/2. These four states are similar in composition 




FIG. 8: Plot of AJ as a function of Xb and x c in the case 
of three mutually interacting electrons. Axis directions are 
reversed from the preceding figure. 




FIG. 9: Plot of A - as a function of Xb and x c in the case 
of three mutually interacting electrons, computed within the 
HM approximation. 



to the four members of the HL 5 = 1/2 subspace. The 
remaining 12 states of the HM 5 — 1/2 subspace consist 
mainly of states with two electrons on a single dot. The 
four 5 = 3/2 states have an energy that is in between 
^i/2,hund an d the energy of the higher-lying 5=1/2 
states. We thus have a situation that is completely anal- 
ogous to the one we encountered in the HL case. It is 
reasonable to project out the 8 low energy states of the 
HM calculation and compare with the HL calculation. 
Figures and show the values of K and J for an 
effective Hamiltonian of the form i|16|) that gives this 8 
dimensional low-energy subspace's spectrum. These fig- 
ures should be compared to Figs. and [7| 

For reasonable parameter values (again, Xf, — 1.0, y = 
1.5), we find in the HL approximation that J = 1.9 for 
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FIG. 10: Plot of J as a function of Xb and x c in the case 
of three mutually interacting electrons, computed within the 
HM approximation. 



two particles (in units of fiu ), J = 1.2 for three parti- 
cles, indicating a change of —36% (or an absolute change 
of —0.7). In the HM approximation, J = 2.6 for two par- 
ticles and J — 1.8 for three particles, indicating a change 
of —31% (or an absolute change of —0.8). Thus, the same 
effect is seen. The absolute value of J is larger in the HM 
case (this is expected since the basis has increased, lead- 
ing to a decrease in the ground state energy E\/2 while 
E3/2 stays constant), but the qualitative HL conclusions 
are well substantiated. 



C. Four Electrons 

The actual calculation for four electron case is more 
involved than that of the three electron case but identical 
in procedure. Details are given in Appendix El The 
resulting quantities K and AK appear as functions of Xb 
and x c in Fig. Illland Fig.^J respectively. Here, AK is 
the value of K minus the value of K obtained by setting 
to zero both three-body (pi, ei) and also four-body (poj 
Pq, eo, e o) matrix elements. 

The behavior of the exchange-interaction constant J 
as a function of Xb and x c (Fig. I13fl is similar to that 
of the three-electron case (Fig. 0). The appearance of 
A J (Fig. El given by subtracting from J the value that 
J would take if the three-body and four-body matrix 
elements were zero) is also reminiscent of A J in the three- 
electron case (Fig. |HJ. On the other hand, J' (Fig. IT5|) 
exhibits rich behavior while A J' fFig. H6|) is qualitatively 
similar in form to A J from the three-electron case. 

The interaction constant J' can be quite significant 
compared to J, which is remarkable and requires atten- 
tion in quantum computer design. In fact, at the point 
Xb = 1, x c = 1.5, our calculation yields J = 0.76 and 
J' = —0.11 so J' I J — —15%, implying substantial four- 
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FIG. 11: Plot of K as a function of xt and x c in the case of 
four mutually interacting electrons. 




FIG. 12: Plot of AK as a function of Xb and x c in the case 
of four mutually interacting electrons. Axis directions are 
reversed from the preceding figure. 

body interactions. We caution, though, that these val- 
ues were obtained within a HL approximation that will 
become inaccurate as Xb decreases and the minima of 
(|26|l get closer together. Our intention is to highlight the 
possible significance of the four-body terms. Such terms 
have been observed experimentally in 3 He, 4Q and CU4O4 
square plaquettes in L^CuOj^i where J' /J was found 
to be ~ 27%. 



V. CONCLUSIONS 

The exchange interaction between localized electrons is 
a basic phenomenon of condensed matter physics, with 
a history that dates back to Heisenberg's pioneering 
worki The details of its behavior are of great signifi- 
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FIG. 13: Plot of J as a function of x b and x c in the case of FIG. 15: Plot of J' as a function of x b and x c in the case of 
four mutually interacting electrons. four mutually interacting electrons. 




FIG. 14: Plot of AJ as a function of x b and x c in the case 
of four mutually interacting electrons. Axis directions are 
reversed from the preceding figure. 




FIG. 16: Plot of AJ' as a function of x b and x c in the case 
of four mutually interacting electrons. Axis directions are 
reversed from the preceding figure. 



cance to quantum information processing using quantum 
dots. Here we have considered the effects that arise when 
three or more electrons, each localized in a low energy or- 
bital on a quantum dot, are simultaneously coupled. We 
have shown that both quantitative and qualitative effects 
arise, due to many-body terms, that modify the stan- 
dard form of the Heisenberg exchange interaction. Most 
significantly, in the case of four coupled electrons, there 
is a four-body interaction that is added to the Heisen- 
berg exchange interaction, and our HL calculations sug- 
gest that it could be strong in physically relevant pa- 
rameter regimes. This possibility needs to be consid- 
ered in electron-spin based quantum computer design, 
because, on the one hand, of the problems it could pro- 
duce when its presence is unwelcome and, on the other 
hand, because of its potential uses in novel designs. In 



other designs as well, the possibility should be considered 
that many-qubit terms could arise in the effective qubit 
Hamiltonian. 
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APPENDIX A: PERTURB ATIVE APPROACH TO 
HEISENBERG EXCHANGE 

Here we summarize the perturbative approach to de- 
riving correcti ons to th e Heisenberg exchange interaction. 
See, e.g., Refs. I7l2fil27l for more details. 

After second quantization of the Coulomb interaction 
Hamiltonian J2J one arrives at the result 

H = ^ £in is + E E *»i°Ia a J* + U E n i^i,-s 

i,s l<j S i 

i<j s,s' 

where a\ s creates an electron with spin s in the ith Wan- 
nier orbital <fi(r— Tj), ni s — a\ s ai S is the number operator, 

hj = J (j)*{r - Ti)h{r)<f>(r - r^dr 

is the hopping energy for i ^ j , Ei = ta is the energy of 
the electron in the ith orbital, 



Wr)IW)l 



l\\2 



■ dvdr' 



is the on-site interaction energy, and 

2 / v (r-ri)0(r-rj)^*(r' - r 3 -)«/>(r' - r<) 

> 



drdr' 



r - r' 



is the customary direct exchange integral (ferromag- 
netic) . 

One now assumes ?7 3> £% , , to ensure that all 
orbitals are singly occupied. One then evaluates H in 
the HL basis 



\si,s 2 , ...,s n ) = a[ Si al 



S2 ns 



a' vac 



in which states differ one from the other only in the dis- 
tribution of spins Si in the orbitals i. The evaluation is 
facilitated by noting that the operators 



are the component of spin 1/2 operators, where a" s , are 
the matrix elements of the Pauli matrices (a = x, y, z ). 
This allows one to rewrite the exchange term as 



E 



.a] s ,a 



j S — 2t/jyS^ ■ Sj — \~ const ^ 



which is the familiar Heisenberg exchange Hamiltonian. 

The contribution of the hopping term to H can be ne- 
glected in the limit tij/U — ► 0. However, when hj/U <C 1 
but nonvanishing, it can be shown, 7 using standard per- 
turbation theory in powers of 1/U, that the effective hop- 
ping Hamiltonian in the HL basis takes the form 



TT _ H h 

~~ U~ ^ IP 



Hi 



El 

U 3 



where Hi, 



s> a is' 



l h = H l <jYjs t ii a is a 3s is tne original hopping 
Hamiltonian, which vanishes in the HL basis. The first 

H 2 

order correction — ^ gives rise to a term of the form 
T! Si<j | ^ | 2 S i • Sj + const, which quantitatively mod- 
ifies (with opposite sign, i.e., antiferromagnctically) the 
Heisenberg Hamiltonian. 6 . However, it is clear that higher 
order terms can contribute multi-spin terms of the form 
we have considered in this paper. It can be shown^& that 
all odd orders vanish, in agreement with our general sym- 

metry argument of Section |H| The term jfi- then gives 
rise to four-spin interactions of the form (Si • Sj)(Sfe • S;), 
proportional to tijtjktkitu/U , with i < j < I, i < k, 
k 7^ j, l. 26 This can be interpreted diagramatically as a 
cycle in which the electrons interchange dots in the order 
i—*l—*k—*j—*i. Thus, perturbation theory shows 
that when t /U 3 is significant, the four-spin interaction 
cannot be neglected. 



APPENDIX B: DETAILS OF CALCULATIONS 
FOR THREE ELECTRONS 

The Hamiltonian @ contains 3 one-body terms h and 
3 Coulomb interaction terms w, and the contribution of 
each term is given in l|27|l - l|29|) . These contributions 
determine the parameters P3, pi, po, £3, ei, and eo that 
appear in equations IfT^ and (fTSft : 



e 3 = (A\h\A) + (B\h\B) + (C\h\C) + (AB\w\AB) + (AC\w\AC) + (BC\w\BC) 
e x = (B\h\A)(A\B) + (A\h\B)(B\A) + (C\h\C) 

+ (BA\w\AB) + (BC\w\AC)(A\B) + (AC\w\BC)(B\A) 
e = (C\h\A)(A\B)(B\C) + (A\h\B)(C\A)(B\C) + (B\h\C)(C\A)(A\B) 

+ (CA\w\AB)(B\C) + (CB\w\AC)(A\B) + (AB\w\BC)(C\A) 

I 



and 



Ps = 1 

pt = (A\B)(B\A) 

po = (C\A)(A\B)(B\C) 



We apply the symmetries of an equilateral triangle, 
(B\A) = (C\A) = (SIC), (B\h\A) = (C\h\A) = (B\h\C), 
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(A\h\A) = (B\h\B) = (C\h\C) to get the matrix elements 
not explicitly listed in (|27[l - (|29[) . All matrix elements 
are functions of hu) , Xb, and x c , so Lq and L\ obtained 
from (|14fl and l|15l) , and K and J obtained from l|17fl and 
((TBj) . are functions of the same quantities. 

APPENDIX C: DETAILS OF CALCULATIONS 
FOR FOUR ELECTRONS 

Here, there is an electron in a Gaussian orbital at each 
of the four potential minima of (|26[) . We use the ana- 

I 



lytical expressions (|27|l - (|29fl to evaluate the many-body 
matrix elements that appear in equations (|20|1 . (|21l) . and 
1|22JI . The Hamiltonian J3J contains four one-body terms 
h and six Coulomb interaction terms w. Taking them all 
into account, we have 



e 4 = (A\h\A) + (B\h\B) + (C\h\C) + (D\h\D) + (AB\w\AB) + (AC\w\AC) 

+ (AD\w\AD) + (BC\w\BC) + (BD\w\BD) + (CD\w\CD), 
e 2 = (B\h\A)(A\B) + (A\h\B)(B\A) + (C\h\C)(B\A)(A\B) + (D\h\D)(B\A)(A\B) 

+ (BA\w\AB) + (BC\w\AC)(A\B) + (BD\w\AD)(A\B) 

+ (AC\w\BC)(B\A) + (AD\w\BD)(B\A) + (CD\w\CD)(B\A)(A\B), 
ex = (A\h\A)(D\B)(B\C)(C\D) +(D\h\B)(B\C)(C\D) + (B\h\C)(D\B)(C\D) 

+ (C\h\D)(D\B)(B\C) + (AD\w\AB)(B\C)(C\D) + (AB\w\AC)(D\B)(C\D) 

+ (AC\w\AD)(D\B)(B\C) + (DB\w\BC)(C\D) + (DC\w\BD)(B\C) 

+ (BC\w\CD)(D\B), 
e = (B\h\A)(A\B)(D\C)(C\D) + (A\h\B)(B\A)(D\C)(C\D) 

+{D\h\C){B\A){A\B){C\D) + (C\h\D)(B\A)(A\B)(D\C) 

+ (BA\w\AB)(D\C)(C\D) + (BD\w\AC)(A\B)(C\D) + (BC\w\AD)(A\B)(D\C) 
+ (AD\w\BC)(B\A)(C\D) + (AC\w\BD)(B\A)(D\C) + (DC\w\CD)(B\A)(A\B), 
^ = (D\h\A)(A\B)(B\C)(C\D)+(A\h\B)(D\A)(B\C)(C\D) 
+ (B\h\C)(D\A)(A\B)(C\D) + (C\h\D)(D\A)(A\B)(B\C) 

+ (DA\w\AB)(B\C)(C\D) + (DB\w\AC)(A\B)(C\D) + (DC\w\AD)(A\B)(B\C) 
+ (AB\w\BC)(D\A)(C\D) + (AC\w\BD)(D\A)(B\C) + (BC\w\CD)(D\A)(A\B). 

The overlap matrix elements are simpler 

Pa = 1 

P2 = (B\A)(A\B) 

Pl = (D\B)(B\C)(C\D) 

po = (B\A)(A\B)(D\C)(C\D) 

p' Q = (D\A)(A\B)(B\C)(C\D). 

Analytical forms are then available for all of the matrix elements of H and all the overlap matrix elements using 
expressions (|27|) - Q29H and using the tetrahedron symmetries 

(B\A) = (C\A) = (D\A) = (B\C) = (B\D) = (C\D) 

(B\h\A) = (C\h\A) = (D\h\A) = (B\h\C) = (B\h\D) = (C\h\D) 

{A\h\A) = (B\h\B) = (C\h\C) = (D\h\D). 



With all of the matrix elements of H and the overlap solving J20J) - (|2"21) and J^SJl - 123 ■ 

matrix elements in hand, we evaluate K 1 J, and J' by 
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